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(54) Adaptive equalization and interpolated timing recovery in a sampled amplitude read channel 
for magnetic recording 

(57) A sampled amplitude read channel for mag- 
netic disk recording which asynchronously samples (24) 
the analog read signal, adaptively equalizes (B103) the 
resulting discrete time sample values according to a tar- 
get partial response, extracts synchronous sample val- 
ues through interpolated timing recovery (B100), and 
detects digital data from the synchronous sample val- 
ues using a Vrterbi sequence detector (34) is disclosed. 
With interpolated timing recovery, the equalizer filter 
and its associated latency are removed from the timing 
recovery loop, thereby allowing a higher order discrete 
time filter and a lower order analog filter. 
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Description 



FIELD OF INVENTION 



s The present invention relates to the control of magnetic disk storage systems for digital computers, particularly to 

a sampled amplitude read channel that employs asynchronous sampling of an analog read signal, adaptive discrete 
time equalization, and interpolated timing recovery. 

CROSS REFERENCE TO RELATED APPLICATIONS AND PATENTS 

10 

This application is related to other co-pending U.S. patent applications, namely application serial numbers 
08/440,515 entitled "Sampled Amplitude Read Channel For Reading User Data and Embedded Servo Data From a 
Magnetic Medium," 08/341,251 entitled "Sampled Amplitude Read Channel Comprising Sample Estimation Equaliza- 
tion, Defect Scanning, Channel Quality, Digital Servo Demodulation, PID Filter for Timing Recovery, and DC Offset 
is Control," 08/313.491 entitled "Improved Timing Recovery For Synchronous Partial Response Recording " and 
08/533,797 entitled "Improved Fault Tolerant Sync Mark Detector For Sampled Amplitude Magnetic Recording " This 
application is also related to several U.S. patents, namely U.S. Pat. No. 5,359,631 entitled "Timing Recovery Circuit for 
Synchronous Waveform Sampling," 5,291,499 entitled "Method and Apparatus for Reduced-Complexity Viterbi-Type 
Sequence Detectors." 5,297,184 entitled "Gain Control Circuit for Synchronous waveform Sampling," 5 329 554 enti- 
re tied "Digital Pulse Detector/ and 5,424,881 entitled "Synchronous Read Channel." All of the above-named patent appli- 
cations and patents are assigned to the same entity, and all are incorporated herein by reference. 

BACKGROUND OF THE INVENTION 
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In magnetic storage systems for computers, digital data serves to modulate the current in a read/write head coil in 
order to write a sequence of corresponding magnetic flux transitions onto the surface of a magnetic medium in concen- 
tric, radially spaced tracks at a predetermined baud rate. When reading this recorded data, the read/write head again 
passes over the magnetic medium and transduces the magnetic transitions into pulses in an analog read signal that 
alternate in polarity. These pulses are then decoded by read channel circuitry to reproduce the digital data. 

Decoding the pulses into a digital sequence can be performed by a simple peak detector in a conventional analog 
read channel or, as in more recent designs, by a discrete time sequence detector in a sampled amplitude read channel 
Discrete time sequence detectors are preferred over simple analog pulse detectors because they compensate for inter- 
symbol interference (ISI) and are less susceptible to channel noise. As a result, discrete time sequence detectors 
increase the capacity and reliability of the storage system. 

There are several well known discrete time sequence detection methods including discrete time pulse detection 
(DPD), partial response (PR) with Viterbi detection, maximum likelihood sequence detection (MLSD), decision-feed- 
back equalization (DFE), enhanced decision-feedback equalization (EDFE), and fixed-delay tree-search with decision- 
feedback (FDTS/DF). 

In conventional peak detection schemes, analog circuitry, responsive to threshold crossing or derivative informa- 
tion, detects peaks in the continuous time analog signal generated by the read head. The analog read signal is "seg- 
mented" into bit cell periods and interpreted during these segments of time. The presence of a peak during the bit cell 
period is detected as a "1" bit, whereas the absence of a peak is detected as a "0" bit. The most common errors in 
detection occur when the bit cells are not correctly aligned with the analog pulse data. Timing recovery, then adjusts 
the bit cell periods so that the peaks occur in the center of the bit cells on average in order to minimize detection errors 
Since timing information is derived only when peaks are detected, the input data stream is normally run length limited 
(RLL) to limit the number of consecutive "0" bits. 

As the pulses are packed closer together on the concentric data tracks in the effort to increase data density detec- 
tion errors can also occur due to intersymbol interference, a distortion in the read signal caused by closely spaced over- 
lapping pulses. This interference can cause a peak to shift out of its bit cell, or its magnitude to decrease resulting in a 
detection error. The ISI effect is reduced by decreasing the data density or by employing an encoding scheme that 
ensures a minimum number of "0" bits occur between "1" bits. For example, a (d,k) run length limited (RLL) code con- 
strains to d the minimum number of "0" bits between "1" bits, and to k the maximum nurrfcer of consecutive "0" bits A 
typical (1 .7) RLL 2/3 rate code encodes 8 bit data words into 12 bit codewords to satisfy the (1 ,7) constraint 

Sampled amplitude detection, such as partial response (PR) with Viterbi detection, allows for increased data den- 
sity by compensating for intersymbol interference and the effect of channel noise. Unlike conventional peak detection 
systems, sampled amplitude recording detects digital data by interpreting, at discrete time instances, the actual value 
of the pulse data. To this end. the read channel comprises a sampling device for sampling the analog read signal and 
a timing recovery circuit for synchronizing the samples to the baud rate (code bit rate). Before sampling the pulses a 
variable gam amplifier adjusts the read signals amplitude to a nominal value, and a tow pass analog filter filters the read 
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signal to attenuate aliasing noise. After sampling, a digital equalizer filter equalizes the sample values according to a 
desired partial response, and a discrete time sequence detector, such as a Viterbi detector, interprets the equalized 
sample values in context to determine a most likely sequence for the digital data (i.e., maximum likelihood sequence 
detection (MLSD)). MLSD takes into account the effect of ISI and channel noise in the detection algorithm, thereby 
5 decreasing the probability of a detection error. This increases the effective signal to noise ratio and, for a given (d.k) 
constraint, allows for significantly higher data density as compared to conventional analog peak detection read chan- 
nels. 

The application of sampled amplitude techniques to digital communication channels is well documented. See Y. 
Kabal and S. Pasupathy, "Partial Response Signaling". IEEE Trans. Commun. Tech., Vol. COM-23, pp.921 -934, Sept. 

io 1975; and Edward A. Lee and David G. Messerschmitt, "Digital Communication", Kluwer Academic Publishers, Boston, 
1990; and G.D. Forney, Jr., "The Viterbi Algorithm", Proc. IEEE, Vol. 61, pp. 268-278, March 1973. 

Applying sampled amplitude techniques to magnetic storage systems is also well documented. See Roy D. Cideci- 
yan, Francois Dolivo, Walter Hirt, and Wolfgang Schott, "APRML System for Digital Magnetic Recording", IEEE Journal 
on selected Areas in Communications, Vol. 10 No. 1, January 1992, pp.38-56; and Wood et al, "Viterbi Detection of 

15 Class IV Partial Response on a Magnetic Recording Channel ", IEEE Trans. Commun., Vol. Com-34, No. 5, pp. 454- 
461, May 1986; and Coker Et al, "Implementation of PRML in a Rigid Disk Drive", IEEE Trans, on Magnetics, Vol. 27, 
No. 6, Nov. 1991; and Carley et al, "Adaptive Continous-Time Equalization Followed By FDTS/DF Sequence Detec- 
tion", Digest of The Magnetic Recording Conference, August 15-17, 1994, pp. C3; and Moon et al, "Constrained-Com- 
plexity Equalizer Design for Fixed Delay Tree search with Decision Feedback", IEEE Trans, on Magnetics, Vol. 30, No. 

20 5, Sept. 1994; and Abbott et al, "Timing Recovery For Adaptive Decision Feedback Equalization of The Magnetic Stor- 
age Channel", G/obecom'90 IEEE Global Telecommunications Conference 1990, San Diego, CA, Nov. 1990, pp. 1794- 
1799; and Abbott et al, "Performance of Digital Magnetic Recording with Equalization and Offtrack Interference", IEEE 
Transactions on Magnetics, Vol. 27, No. 1 , Jan. 1991 ; and Cioffi et al, "Adaptive Equalization in Magnetic-Disk Storage 
Channels", IEEE Communication Magazine, Feb. 1990; and Roger Wood, "Enhanced Decision Feedback Equaliza- 

25 tion". lntermag'90. 

Similar to conventional peak detection systems, sampled amplitude detection requires timing recovery in order to 
correctly extract the digital sequence. Rather than process the continuous signal to align peaks to the center of bit cell 
periods as in peak detection systems, sampled amplitude systems synchronize the pulse samples to the baud rate. In 
conventional sampled amplitude read channels, timing recovery synchronizes a sampling clock by minimizing an error 
30 between the signal sample values and estimated sample values. A pulse detector or slicer determines the estimated 
sample values from the read signal samples. Even in the presence of ISI the sample values can be estimated and, 
together with the signal sample values, used to synchronize the sampling of the analog pulses in a decision-directed 
feedback system. 

A phase-locked-loop (PLL) normally implements the timing recovery decision-directed feedback system. The PLL 
35 comprises a phase detector for generating a phase error based on the difference between the estimated samples and 
the read signal samples. A PLL loop filter filters the phase error, and the filtered phase error operates to synchronize 
the channel samples to the baud rate. Conventionally, the phase error adjusts the frequency of a sampling clock which 
is typically the output of a variable frequency oscillator (VFO). The output of the VFO controls a sampling device, such 
as an analog-to-digital (A/D) converter, to synchronize the sampling to the baud rate. 
40 As mentioned above, sampled amplitude read channels also commonly employ a discrete time equalizer filter to 
equalize the sample values into a desired partial response (PR4, EPR4, EEPR4, etc.) before sequence detection. To 
this end, adaptive algorithms have been applied to compensate in real time for parameter variations in the recording 
system and across the disk radius. For example, U.S. Pat. No. 5,381,359 entitled "Adaptation and Training of Digital 
Finite Impulse Response Filter Within PRML Sampling Data Detection Channel", discloses an adaptive equalizer filter 
45 that operates according to a well known least mean square (LMS) algorithm, 

W k+1 = W k - M .e k .X k , 

where represents a vector of filter coefficients; n is a programmable gain; e k represents a sample error between the 
so filter's actual output and a desired output; and represents a vector of sample values from the filter input. In other 
words, the LMS adaptive equalizer filter is a closed loop feedback system that attempts to minimize the mean squared 
error between an actual output of the filter and a desired output by continuously adjusting the filters coefficients to 
achieve an optimum frequency response. 

A problem associated with adaptive equalizer filters in sampled amplitude read channels is that the timing recovery 
55 and gain control loops can interfere with the adaptive feedback loop, thereby preventing the adaptive equalizer filter 
from converging to an optimal state. This non-convergence is manifested by the filter's phase and gain response drifting 
as it competes with the timing and gain control loops. An article by J.D. Coker et al. entitled "Implementation of PRML 
in a Rigid Disk Drive", published in IEEE Transactions on Magnetics, vol. 27, No. 6, Nov. 1991, suggests a three tap 
transversal filter comprising a fixed center tap and symmetric side taps in order to constrain the phase response of the 
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equalizer filter except in terms of a fixed group delay. Constraining the phase response of the adaptive equalizer in this 
manner, however, is a very sub-optimal method for attenuating interference from the timing recovery and gain control 
loops. Furthermore, it significantly reduces control over the adaptive filter's phase response, thereby placing the burden 
of phase compensation on the analog equalizer. 

Yet another problem associated with conventional adaptive equalizer filters is an inherent limitation on its order (i.e., 
the number of coefficients): because the adaptive equalizer is inside the timing recovery feedback loop, its order must 
be limited to minimize the associated transport delay. Compensating for the def iciencies of the discrete time equalizer 
requires filtering the analog read signal with a higher order analog equalizer prior to the timing recovery loop, which is 
undesirable. 

_ There is, therefore, a need for an adaptive, discrete time equalizer filter in a sampled amplitude read channel hav- 
ing an improved method for constraining the phase and gain response in order to minimize interference from the timing 
recovery and gain control loops. A further aspect of the present invention is to remove the adaptive equalizer, and its 
associated latency, from the timing recovery loop, thereby allowing a higher order discrete time filter and a simplified 
analog filter. 

SUMMARY OF THE INVENTION 

A sampled amplitude read channel for magnetic disk recording which asynchronously samples the analog read sig- 
nal, adaptively equalizes the resulting discrete time sample values according to a target partial response, extracts syn- 
chronous sample values through interpolated timing recovery, and detects digital data from the synchronous sample 
values using a Viterbi sequence detector is disclosed. To minimize interference from the timing and gain control loops, 
the phase and magnitude response of the adaptive equalizer filter are constrained at a predetermined frequency using 
an optimal orthogonal projection operation as a modification to a least mean square (LMS) adaptation algorithm. Fur- 
ther, with interpolated timing recovery, the equalizer filter and its associated latency are removed from the timing recov- 
ery loop, thereby allowing a higher order discrete time filter and a lower order analog filter. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The above and other aspects and advantages of the present invention will be better understood by reading the fol- 
lowing detailed description of the invention in conjunction with the drawings, wherein: 

FIG. 1 is a block diagram of a conventional sampled amplitude recording channel. 

FIG. 2A shows an exemplary data format of a magnetic disk having a plurality of concentric tracks comprised of a 
plurality of user data sectors and embedded servo data sectors. 
FIG. 2B shows an exemplary format of a user data sector. 

FIG. 3 is a block diagram of the improved sampled amplitude read channel of the present invention comprising 
interpolated timing recovery for generating interpolated sample values and a synchronous data clock for clocking 
operation of a discrete time sequence detector and peak detector. 

FIG. 4A is a detailed block diagram of the prior art sampling timing recovery comprising a sampling VFO. 
FIG. 4B is a detailed block diagram of the interpolating timing recovery of the present invention comprising an inter- 
polator. 

FIG. 5 illustrates the channels samples in relation to the interpolated baud rate samples for the acquisition pream- 
ble. 

FIG. 6 shows an FIR fitter implementation for the timing recovery interpolator. 
FIG. 7 depicts a cost reduced implementation for the timing recovery interpolator. 

FIG. 8A is a block diagram of a conventional adaptive, discrete time equalizer filter in a sampled amplitude read 
channel. 

FIG. 8B shows the adaptive, discrete time equalizer of the present invention. 

FIG. 8C shows an alternative embodiment for the adaptive, discrete time equalizer of the present invention. 

FIG. 9A illustrates the present invention adaptive filter's gain response constrained at a normalized frequency of 

1/4T 

FIG. 9B shows the present invention adaptive filter's phase response constrained at a normalized frequency of 
1/4T 

FIG. 1 0 illustrates operation of an orthogonal projection operation of the present invention for constraining the gain 
and phase response of the adaptive filter. 

FIG. 1 1 A shows an implementation for a reduced cost orthogonal projection operation. 
FIG. 1 1 B shows an alternative embodiment for the reduced cost orthogonal projection operation. 
FIG. 1 1C illustrates an implementation for a gradient averaging circuit used in the reduced cost orthogonal projec- 
tion operation of FIG. 1 1 B. 
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FIG. 11D shows yet another alternative embodiment for the reduced cost orthogonal projection operation of the 
present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

5 

Conventional Sampled Amplitude Read Channel 

Referring now to FIG. 1 , shown is a detailed block diagram of a conventional sampled amplitude read channel. Dur- 
ing a write operation, either user data 2 or preamble data from a data generator 4 (for example 2T preamble data) is 

io written onto the media. An RLL encoder 6 encodes the user data 2 into a binary sequence b(n) 8 according to an RLL 
constraint. A precoder 10 precodes the binary sequence b(n) 8 in order to compensate for the transfer function of the 
recording channel 18 and equalizer filters to form a precoded sequence ~b(n) 12. The precoded sequence ~b(n) 12 is 
converted into symbols a(n) 16 by translating 14 ~b(N) = 0 into a(N) = -1, and ~b(N) = 1 into a(N) = +1 . Write circuitry 
9, responsive to the symbols a(n) 1 6, modulates the current in the recording head coil at the baud rate IfT to record the 

is binary sequence onto the media. A frequency synthesizer 52 provides a baud rate write clock 54 to the write circuitry 9 
and is adjusted by a channel data rate signal (CDR) 30 according to the zone the recording head is over. 

When reading the recorded binary sequence from the media, timing recovery 28 first locks to the write frequency 
by selecting, as the input to the read channel, the write clock 54 through a multiplexor 60. Once locked to the write fre- 
quency, the multiplexor 60 selects the signal 1 9 from the read head as the input to the read channel in order to acquire 

20 an acquisition preamble recorded on the disk prior to the recorded user data. A variable gain amplifier 22 adjusts the 
amplitude of the analog read signal 58, and an analog filter 20 provides initial equalization toward the desired response 
as well as attenuating aliasing noise. A sampling device 24 samples the analog read signal 62 from the analog filter 20, 
and a discrete time equalizer filter 26 provides further equalization of the sample values 25 toward the desired 
response. In partial response recording, for example, the desired response is often selected from Table 1 . 

25 The discrete equalizer filter 26 may be implemented as a real-time adaptive filter which compensates for parameter 

variations over the disk radius (i.e.. zones), disk angle, and environmental conditions such as temperature drift. To this 
end. the filter 26 receives estimated sample values B143 generated by the timing recovery circuit 28; the estimated 
samples being input into an adaptive feedback loop and used to generate sample errors. The adaptive feedback loop 
conventionally employs a least mean square (LMS) algorithm to adapt the filter coefficients (i.e. , it adapts the frequency 

30 and phase response of the filter) until a minimum sample error is achieved. Operation of a conventional adaptive equal- 
izer filter is discussed in greater detail below. 

After equalization, the equalized sample values 32 are applied to a decision directed gain control 50 and timing 
recovery 28 circuit for adjusting the amplitude of the read signal 58 and the frequency and phase of the sampling device 
24, respectively. Timing recovery adjusts the frequency of sampling device 24 over line 23 in order to synchronize the 

35 equalized samples 32 to the baud rate. Frequency synthesizer 52 provides a course center frequency setting to the tim- 
ing recovery circuit 28 over line 64 in order to center the timing recovery frequency over temperature, voltage, and proc- 
ess variations. The channel data rate (CDR) 30 signal adjusts a frequency range of the synthesizer 52 according to the 
data rate for the current zone. Gain control 50 adjusts the gain of variable gain amplifier 22 over line 21 in order to match 
the magnitude of the channel s frequency response to the desired partial response. 

40 The equalized samples Y(n) 32 are also sent to a discrete time sequence detector 34, such as a maximum likeli- 
hood (ML) Viterbi sequence detector, which detects an estimated binary sequence A b(n) 33 from the sample values. An 
RLL decoder 36 decodes the estimated binary sequence A b(n) 33 from the sequence detector 34 into estimated user 
data 37. A data sync detector 66 detects the sync mark 70 (shown in FIG. 2B) in the data sector 15 in order to frame 
operation of the RLL decoder 36. In the absence of errors, the estimated binary sequence A b(n) 33 matches the 

45 recorded binary sequence b(n) 8, and the decoded user data 37 matches the recorded user data 2. 

Data Format 

FIG. 2A shows an exemplary data format of a magnetic media comprising a series of concentric data tracks 13 
bo wherein each data track 1 3 comprises a plurality of sectors 1 5 with embedded servo wedges 17. A servo controller (not 
shown) processes the servo data in the servo wedges 17 and, in response thereto, positions the read/write head over 
a desired track. Additionally, the servo controller processes servo bursts within the servo wedges 17 to keep the head 
aligned over a centerline of the desired track while writing and reading data. The servo wedges 17 may be detected by 
a simple discrete time pulse detector or by the discrete time sequence detector 34. If the sequence detector 34 detects 
55 the servo data, then the format of the servo wedges 17 includes a preamble and a sync mark, similar to the user data 
sectors 15. 

FIG. 2B shows the format of a user data sector 15 comprising an acquisition preamble 68, a sync mark 70, and 
user data 72. Timing recovery uses the acquisition preamble 68 to acquire the correct sampling frequency and phase 
before reading the user data 72, and the sync mark 70 demarks the beginning of the user data 72 (see co-pending U.S. 
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patent application 08/313,491 entitled "Improved Timing Recovery For Synchronous Partial Response Recording"). 

To increase the overall storage density, the disk is partitioned into an outer zone 11 comprising fourteen data sec- 
tors per track, and an inner zone 27 comprising seven data sectors per track. In practice, the disk is actually partitioned 
into several zones with a different number of sectors in each zone, and the data recorded and detected at a different 
data rate in each zone. 

Improved Sampled Amplitude Read Channel 

FIG. 3 shows the improved sampled amplitude read channel of the present invention wherein the conventional 
sampled timing recovery 28 of FIG. 1 has been replaced by interpolated timing recovery B100. In addition the write 
frequency synthesizer 52 generates a baud rate write clock 54 applied to the write circuitry 9, and an asynchronous 
read clock 54 for clocking the sampling device 24. the discrete time equalizer filter 26, and the interpolated timing recov- 
ery B1 00 at a frequency relative to the current zone (CDR 30). In an alternative embodiment, a first frequency synthe- 
sizer generates the write clock, and a second frequency synthesizer generates the read clock. 

The discrete equalizer filter B103 is real-time adaptive, receiving interpolated sample values B102 and estimated 
sample values B143 from the interpolated timing recovery circuit B100 for use in a modified least mean square (LMS) 
algorithm wh.ch constrains the filter's gain and phase response according to the present invention, the details of which 
are set forth below. 

The interpolated timing recovery B100 interpolates the equalized sample values 32 to generate interpolated sam- 
ple values B102 substantially synchronized to the data rate of the current zone. A discrete time sequence detector 34 
detects an estimated binary sequence 33 representing the user data from the interpolated sample values B102 The 
interpolated timing recovery B100 circuit generates a synchronous data clock B104 for clocking operation of the gain 
control 50, discrete time sequence detector 34, sync mark detector 66 and RLL decoder 36. 

25 Conventional Timing Recovery 

An overview of the conventional sampling timing recovery 28 of FIG. 1 is shown in FIG. 4A. The output 23 of a var- 
iable frequency oscillator (VFO) B164 controls the sampling clock of a sampling device 24 which is typically an analog- 
to-digrtal converter (A/D) in digital read channels. A multiplexor B159 selects the unequalized sample values 25 during 
so acquisition and the equalized sample values 32 during tracking, thereby removing the discrete equalizer filter 26 from 
the timing loop during acquisition in order to avoid its associated latency. A phase error detector B155 generates a 
phase error in response to the sample values received over line B149 and estimated sample values ~Y k from a sample 
value estimator B141. such as a slicer in a d=0 PR4 read channel, over line B143. A loop filter B160 filters the phase 
error to generate a frequency offset Af B167 that settles to a value proportional to a frequency difference between the 
samphng clock 23 and the baud rate. The frequency offset Af B167. together with the center frequency control signal 
64 from the frequency synthesizer 52, adjust the sampling clock 23 at the output of the VFO B164 in order to synchro- 
nize the sampling to the baud rate. 

A zero phase start B162 circuit suspends operation of the VFO B164 at the beginning of acquisition in order to min- 
imize the initial phase error between the sampling clock 23 and the read signal 62. This is achieved by disabling the 
VFO B164. detecting a zero crossing in the analog read signal 62. and re-enabling the VF0 164 after a predetermined 
delay between the detected zero crossing and the first baud rate sample. 

IT* H? mated samp,e values 8143 at The o^P" 1 of th e slicer B141 are also input into the discrete time equalizer 
filter 26 of FIG. 1 for use in a conventional least mean square (LMS) adaptation algorithm as is described in more detail 
below. 
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Interpolated Timing Recovery 

The interpolated timing recovery B100 of the present invention is shown in FIG. 4B. The VFO B164 in the conven- 
tional tiirung recovery of FIG. 4A is replaced with a modulo-Ts accumulator B120 and an interpolator B122. In addition 
an expected sample value generator B151. responsive to interpolated sample values B102. generates expected sam- 
ples Y^ used by the phase error detector B155 to compute the phase error during acquisition A multiplexor B153 
selects the estimated sample values ~Y k+t from the slicer B141 for use by the phase error detector B155 during track- 
ing. The data clock B104 is generated at the output of an AND gate B126 in response to the sampling clock 54 and a 
mask signal B1 24 from the modulo-Ts accumulator B1 20 as discussed in further detail below. The phase error detector 
B155 and the sheer B141 process interpolated sample values B102 at the output of the interpolator B122 rather than 
the channel sample values 32 at the output of the discrete equalizer filter 26 as in FIG. 4A. A PID loop filter B161 con- 
trols the closed loop frequency response similar to the loop filter B160 of FIG. 4A. 

The interpolated sample values Y k+t B102 and the estimated sample values ~Y k „ from the slicer B141 are input 
into the adaptive, discrete equalizer filter B103 of FIG. 3 for use by a modified least mean square (LMS) algorithm, the 
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details of which are set forth below. 

in the interpolated timing recovery of the present invention, locking a VFO to a reference frequency before acquir- 
ing the preamble is no longer necessary; multiplexing 60 the write clock 54 into the analog receive filter 20 (as in FIG. 
1) is not necessary. Further, the sampling device 24 and the discrete equalizer filter 26, together with their associated 

5 delays, have been removed from the timing recovery loop; it is not necessary to multiplex B159 around the equalizer 
filter 26 between acquisition and tracking. However, it is still necessary to acquire a preamble 68 before tracking the 
user data 72. To this end, a zero phase start circuit B163 minimizes the initial phase error between the interpolated 
sample values and the baud rate at the beginning of acquisition similar to the zero phase start circuit B162 of FIG. 4A. 
However, rather than suspend operation of a sampling VFO B1 64, the zero phase start circuit B1 63 for interpolated tim- 

w ing recovery computes an initial phase error x from the A/D 24 sample values 25 and loads this initial phase error into 
the modulo-Ts accumulator B120. 

For more details concerning the PID loop filter B1 61 , phase error detector B1 55, expected sample generator B1 51 , 
and slicer B141 . refer to the above referenced co-pending U.S. patent applications "Sampled Amplitude Read Channel 
Comprising Sample Estimation Equalization, Defect Scanning, Channel Quality, Digital Servo Demodulation, PID Filter 

75 for Timing Recovery, and DC Offset Control" and "Improved Timing Recovery For Synchronous Partial Response 
Recording " A detailed description of the modulo-Ts accumulator B120, data clock B104, and interpolator B122 is pro- 
vided in the following discussion. 

Interpolator 

20 

The interpolator B1 22 of FIG. 4B is understood with reference to FIG. 5 which shows a sampled 2T acquisition pre- 
amble signal B200 The target synchronous sample values B102 are shown as black circles and the asynchronous 
channel sampte values 32 as vertical arrows. Beneath the sampled preamble signal is a timing diagram depicting the 
corresponding timing signals for the sampling clock 54, the data clock B104 and the mask signal B124. As can be seen 
25 in FIG 5. the preamble signal B200 is sampled slightly faster than the baud rate (the rate of the target values). 

The function of the interpolator is to estimate the target sample value by interpolating the channel sample values. 
For illustrative purposes, consider a simple estimation algorithm, linear interpolation: 

Y(N-1) = x(N-1) + t * (x(N) - x(N-1)); where: (1) 

30 

x(l\M) and *(N) are the channel samples surrounding the target sample; and t is an interpolation interval proportional 
to a time dittei ence between the channel sample value x(N-1) and the target sample value. The interpolation interval t 
is generated at the output of modulo-Ts accumulator B120 which accumulates the frequency offset signal Af B167 at 
the output o< the PID loop filter B161 : 

35 

x = (£A0 MOD TS; where: (2) 

Ts is the sampling peuod of the sampling clock 54. Since the sampling clock 54 samples the analog read signal 62 
slightly tastef rvin the baud rate, it is necessary to mask the data clock every time the accumulated frequency offset Af , 

40 integer drvded by Ts. increments by 1 . Operation of the data clock B104 and the mask signal B124 generated by the 
modulo-Ts accumulator B120 is understood with reference to the timing diagram of FIG. 5. 

Assuming me interpolator implements the simple linear equation (1) above, then channel sample values B202 and 
B204 are usee to gt?oei ate the interpolated sample value corresponding to target sample value B206. The interpolation 
interval i B208 is generated according to equation (2) above. The next interpolated sample value corresponding to the 

45 next target value B210 is computed from channel sample values B204 and B212. This process continues until the inter- 
polation mtef vai \ B2 1 4 would be greater than Ts except that it "wraps" around and is actually x B21 6 (i.e. , the accumu- 
lated frequency onset \'. integer divided by Ts, increments by 1 causing the mask signal B1 24 to activate). At this point, 
the data etoefc B104 is masked by mask signal B124 so that the interpolated sample value corresponding to the target 
sample value B220 is computed from channel sample values B222 and B224 rather than channel sample values B218 

so and B222 

The simple linear interpolation of equation (1) will only work if the analog read signal is sampled at a much higher 
frequency than the baud rate. This is not desirable since operating the channel at higher frequencies increases its com- 
plexity and cost Therefore, in the preferred embodiment the interpolator B122 is implemented as a filter responsive to 
more than two channel samples to compute the interpolated sample value. 
55 The ideal discrete time phase interpolation filter has a flat magnitude response and a constant group delay of x: 

C T (e jc W^ (3) 

which has an ideal impulse response: 
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sinc(7c-(n-x/T s )). (4 ) 

Unfortunately, the above non-causal infinite impulse response (4) cannot be realized. Therefore, the impulse response 
of the interpolation filter is designed to be a best fit approximation of the ideal impulse response (4). This can be accom- 
plished by minimizing a mean squared error between the frequency response of the actual interpolation filter and the 
frequency response of the ideal interpolation filter (3). This approximation can be improved by taking into account the 
spectrum of the input signal, that is, by minimizing the mean squared error between the input spectrum multiplied by 
the actual interpolation spectrum and the input spectrum multiplied by the ideal interpolation spectrum: 

C T (e iw )X(e iG) ) - C T (e ja) )X(e jo ); where: (5 ) 

C^eh is the spectrum of the actual interpolation filter; and Xfe^) is the spectrum of the input signal. From equation 
(5), the mean squared error is represented by:_^ 



E ' = ^J |C * (eK °>-^ (6) 



20 X(e^) is the spectrum of the read channel (e.g., PR4, EPR4, EEPR4 of Table 1 or some other partial response spec- 
trum). 

In practice, the above mean squared error equation (6) is modified by specifying that the spectrum of the input sig- 
nal is bandlimitted to some predetermined constant 0 < <*> < an where 0 < a < 1 ; that is: 

25 |X(e ico )| = 0,for |<o|> an. 

Then equation (6) can be expressed as: 



30 



Ej.« = ^ J lC t (e j ") - en 2 \X(ehl 2 d^ (7 ) 



The solution to the minimization problem of equation (7) involves expressing the actual interpolation filter in terms of its 
35 coefficients and then solving for the coefficients that minimize the error in a classical mean-square sense. 
The actual interpolation filter can be expressed as the FIR polynomial: 



n=R-1 

ico, 



40 



C t (e ,ffl )= £ C T (n)e " >n ; where: {8 ) 



n=-R 



45 



50 



2R is the number of taps in each interpolation filter and the sample period Ts has been normalized to 1. A mathematical 
derivation for an interpolation filter having an even number of coefficients is provided below. It is within the ability of 
those skilled in the art to modify the mathematics to derive an interpolation filter having an odd number of coefficients. 
Substituting equation (8) into equation (7) leads to the desired expression in terms of the coefficients C T (n): 

tin n=R-i 

E *.« =2^ J I £ C t (n)e jmn - e^Wtfda,. (9) 

-an n=-R 

The next step is to take the derivatives of equation (9) with respect to the coefficients C T (n) and set them to zero: 

55 2 

^fo = 0forn o-R 0.1 R-1. (10) 

After careful manipulation, equation (10) leads to: 
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J \[ Z C T (n)cos(a>(n 0 - n))l - cos( «(n D + T))l|X(e jQ> )| 2 d<i> = 0 



(11) 



for n 0 = -R. .... 0, 1, R-1 
Defining 4>(r) as: 



10 



<K0= J |X(e jto )| 2 cos(<or) da> 



(12) 



15 



and substituting equation (12) into equation (11) gives: 



n=R-1 

X C t (n)<Kn-n 0 ) = (Kn 0 + T) 



(13) 



20 



25 



30 



forn Q = -R 0, 1 R-1 . 

Equation (13) defines a set of 2R linear equations in terms of the coefficients C T (n). Equation (13) can be expressed 
more compactly in matrix form: 



C T is a column vector of the form: 



0> T is a Toeplitz matrix of the form: 



d> T C T = <b x ; where: 



C T = [c T (-R) c T (0), ....c^R-l)] 1 



35 



40 



Ot = 



♦ (0) 

<Mi> 

<t>(2R-l) 



<J>< 1 > 
<M0) 



<J>(2R-1 



♦ (0) 



45 



SO 



and <I> T is a column vector of the form: 



* T = [*(-R + t) <K*), 4>0 + x) <KFM + x)] \ 



The solution to equation (14) is: 



C T = 4> T " <D T ; where: 



<t> T ~^ is an inverse matrix that can be solved using well known methods. 



(14) 



(15) 



Table B2 shows example coefficients C x (n) calculated from equation (15) with 2R=6, a=0.8 and X(e ico )=PR4 . The 
implementation of the six tap FIR filter is shown in FIG. 6. A shift register B250 receives the channel samples 32 at the 
sampling clock rate 54. The filter coefficients C T (n) are stored in a coefficient register file B252 and applied to corre- 
55 sponding multipliers according to the current value of x B128. The coefficients are multiplied by the channel samples 
32 stored in the shift register B250. The resulting products are summed B254 and the sum stored in a delay register 
B256. The coefficient register file B252 and the delay register B256 are clocked by the data clock B104 to implement 
the masking function described above. 

In an alternative embodiment not shown, a plurality of static FIR filters, having coefficients that correspond to the 
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different values of x, filter the sample values in the shift register B250. Each filter outputs an interpolation value, and the 
current value of the interpolation interval x B1 28 selects the output of the corresponding filter as the output B1 02 of the 
interpolator B122. Since the coefficients of one filter are not constantly updated as in FIG. 6, this multiple filter embod- 
iment increases the speed of the interpolator B122 and the overall throughput of the read channel. 

Cost Reduced Interpolator 

Rather than store all of the coefficients of the interpolation filters in memory, in a more efficient, cost reduced imple- 
mentation the coefficient register file B252 of FIG. 6 computes the filter coefficients C t (n) in real time as a function of x. 
For example, the filter coefficients C t (n) can be computed in real time according to a predetermined polynomial in x 
(see, for example. U.S. Pat. No. 4,866,647 issued to Farrow entitled, "A Continuously Variable Digital Delay Circuit," the 
disclosure of which is hereby incorporated by reference). An alternative, preferred embodiment for computing the filter 
coefficients in real time estimates the filter coefficients according to a reduced rank matrix representation of the coeffi- 
cients. 

15 The bank of filter coefficients stored in the coefficient register file B252 can be represented as an MxN matrix Am xN . 
where N is the depth of the interpolation filter (i.e.. the number of coefficients C T (n) in the impulse response computed 
according to equation (15)) and M is the number of interpolation intervals (i.e., the number of x intervals). Rather than 
store the entire A„ kN matrix in memory, a more efficient, cost reduced implementation is attained through factorization 
and singular value decomposition (SVD) of the A^ xN matrix. 

20 Consider that the matrix can be factored into an F MxN and matrix, 

^ MXN = ^ MxN * G NxN . 

Then a reduced rank approximation of the A MxN matrix can be formed by reducing the size of the F MxN and G NxN matri- 
x's ces by replacing N with L where L<N and, preferably. L«N. Stated differently, find the F MxL and matrices whose 
product best approximates the A MxN matrix. 

A MxN ~ F Mx| _ • G UN . 

30 The convolution pi ocess of the interpolation filter can then be carried out, as shown in FIG. 7, by implementing the G LxN 
matrix as a bank of FIR fitters B260 connected to receive the channel sample values 32, and the F MxL matrix imple- 
mented as a looM) table B262 indexed by x B128 (as will become more apparent in the following discussion). Those 
skilled m the art will iecogrnze that, in an alternative embodiment, the A MxN matrix can be factored into more than two 
matrices (i e . A * FGH ) 



35 



40 



50 



The preferred method tor finding the F MxL and G^m matrices is to minimize the following sum of squared errors: 



M N 



Z Z < A jn ' < F MxL # G|_xN)jn) 2 (16) 



J = 



The solution to equaton ( 16) can be derived through a singular value decomposition of the A^ matrix, comprising the 

steps of . 

45 1 performing an SVD on the A MxN matrix which gives the following unique factorization (assuming M>N): 

A mxn = u MxN * D NxN * v nxn where: 



U MaN *aMiN unitary matrix; 

D N>N is a N*N diagonal matrix {a^ <y 2 , .... a N } where <7j are the singular values of A^, and <n > <7o.„ > a w > 
0;and * N 

V NxN ts a NxN unrtary matrix; 

2. selecting a predetermined L number of the largest singular values a to generate a reduced size diagonal matrix 

& D LxL : 
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w 



D LxL =Diag {a,,a 2 , a L } = 



Oj 0 ... 
0 a 2 0... 



0 



0 



. 0 

0 cTt 



25 



3. extracting the first L columns from the U MxN matrix to form a reduced U MxL matrix: 



U, 



MxL ~ 




... u 



1,N 



u 



M,N 



35 



40 



4. extracting the first L rows from the V NxN matrix to form a reduced V Lx n matrix: 



V, 



LxN - 




5. defining the F MxL and G|_ x n matrices such that: 

F MxL * G UN " U MxL " D LxL * V LxN ~ A MxN 

so (for example, let F MxL = U MxL • D ^ and G UN = V ^ ). 

In the above cost reduced polynomial and reduced rank matrix embodiments, the interpolation filter coefficients 
C t (n) are computed in real time as a function of x; that is. the filter's impulse response h(n) is approximated according to: 



is 



h(n.x) = c T (n) = ^Gj(n) • f(i,x): where: 



(17) 
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f(i,x) is a predetermined function in x (e.g., polynomial in t or t indexes the above F MxL matrix); L is a degree which 
determines the accuracy of the approximation (e.g., the order of the polynomial or the column size of the above F MxL 
matrix); and Gj(n) is a predetermined matrix (e.g.. the coefficients of the polynomial or the above G UN matrix). As L 
increases, the approximated filter coefficients C t (n) of equation (17) tend toward the ideal coefficients derived from 
5 equation (1 5). It follows from equation (1 7) that the output of the interpolation filter Y(x) can be represented as: 

N L 

Y(x)= £ u (™) Z G i< n ) ' HI*) (18) 

n= 1 i = 1 

10 

where U(x) are the channel sample values 32 and N is the number of interpolation filter coefficients C T (n). 

Referring again to FIG. 6, the coefficient register file can compute the interpolation filter coefficients C x (n) according 
to equation (1 7) and then convolve the coefficients C x (n) with the channel samples U(x) 32 to generate the interpolated 
sample values B102 synchronized to the baud rate. However, a more efficient implementation of the interpolation fitter 
15 can be achieved by rearranging equation (18): 

l N 

Y(x) = £f(i,x) £ Grfn) • U(x-n) (19) 

i = 1 n= 1 

20 

FIG. 7 shows the preferred embodiment of the interpolation filter according to equation (19). In the polynomial 
embodiment, the function of x is a polynomial in t, and the matrix Gj(n) are the coefficients of the polynomial. And in the 
reduced rank matrix embodiment, the function of t is to index the above F MxL matrix B262, and the second summation 
25 in equation (19), 



£ G i(n) * U(x-n) 

n= 1 

30 

is implemented as a bank of FIR filters B260 as shown in FIG. 7. Again, in equation (19) L is the depth of the approxi- 
mation function f(i, T ) (e.g., the order of the polynomial or the column size of the above F MxL matrix) and N is the depth 
of the interpolation filter's impulse response (i.e., the number of coefficients in the impulse response). It has been deter- 
35 mined that N=8 and U3 provides the best performance/cost balance; however, these values may increase as IC tech- 
nology progresses and the cost per gate decreases. 

Conventional Adaptive Equalizer 

40 FIG. 8A illustrates a prior art adaptive, discrete time equalizer that operates according to the well known least mean 
square (LMS) algorithm, 

=W k -n-e k -X k , 

45 or alternatively, 

W k+ i = W k - M -X k .e k 

where represents a vector of FIR filter coefficients; n is a programmable gain; represents a sample error (or vec- 
so tor of sample errors ej between the FIR filter's actual output and a desired output; and X„ represents a vector of sam- 
ples values (or a scalar from the FIR filter input To better understand operation of the present invention, the second 
representation of the LMS algorithm is used throughout this disclosure. 

The desired filter output is the estimated sample values ~ Y k at the output of slicer B141 . The estimated sample val- 
ues ~Y k are subtracted from the FIR filter's output Y k to generate the sample error e^ The LMS algorithm attempts to 
55 minimize the sample error in a least mean square sense by adapting the FIR filter coefficients; that is. it adjusts the FIR 
filter's gain and phase response so that the overall channel response adequately tracks the desired partial response 
(e.g., PR4, EPR4, EEPR4, etc.). 

As previously mentioned, interference from the timing recovery 28 and gain control 50 loops can prevent the adap- 
tive, discrete time equalizer 26 from converging to an optimal state. For example, a phase adjustment in the adaptive 
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filter 26 can affect the sampling phase error for timing recovery 28. Timing recovery 28 compensates for the filter's 
phase adjustment by adjusting its sampling phase; this adjustment can result in yet another phase adjustment by the 
adaptive equalizer 26. Thus, the phase response of the adaptive equalizer may never converge. Similarly the gain con- 
trol loop 50 can interfere with the gain response of the adaptive filter 26 and prevent it from converging. 

5 

Constrained Adaptive Equalizer 

FIG. 8B illustrates operation of the adaptive, discrete time equalizer filter B1 03 of the present invention. Sample val- 
ues from the A/D 24 are input over line 25 to a discrete time FIR filter C1 00 comprising a predetermined number of coef- 

io_ ficients, the values of which determine the filter's gain and phase response. Because the FIR filter C100 operates on 
the sample values prior to the interpolated timing recovery loop B1 00, its order can be increased over the prior art with- 
out adversely affecting the latency of timing recovery (i.e., the number of filter coefficients can be increased). 

The output Y k 32 of the FIR filter C100 is input into the interpolator B122 for generating the interpolated sample 
values Y k+T B102. The interpolated sample values Y k+t B102 are input into a slicer B141 (FIG. 4B) which generates 

is estimated sample values ~Y k+T . The estimated sample values — Y k+T are subtracted from the interpolated sample val- 
ues Y k+T at adder C102 to generate a sample error value e k+T C104 that is synchronized to the baud rate rather than 
the sample rate. Because the LMS algorithm operates on sample values X k at the sample rate, it is necessary to con- 
vert the error value e^ C104 into an error value e k C112 synchronous to the sample rate. This is accomplished by an 
interpolation circuit C106 which computes an interpolated error value eK C112 from the baud rate error values e K+T 

20 CI 04. Preferably, the error value interpolation circuit C1 06 is implemented as a first order linear interpolation, but it may 
be a simple zero order hold, or a more complex interpolation filter as described above. 

FIG. 8C shows an alternative embodiment for generating the error value e^ As illustrated, the estimated ideal sam- 
ple values ~Y k+T from the slicer B141 are interpolated by the interpolation circuit C106 to generate estimated ideal sam- 
ple values ~Y k C120 which are subtracted C102 from the equalized sample values Y k 32 at the output of the FIR filter 

25 C1 00 to generate the error value e k 

In both embodiments, the error value e k C112 is input into a modified LMS circuit C1 14 which computes updated 
filter coefficients C1 16 according to. 

Wk*i = Wk-M-Pv lV2 x .(X k -e k ) 

30 

where the operation Pv 1 v 2 - L is an orthogonal projection operation which constrains the gain and phase response of the 
FIR filter C1 00 in order to attenuate interference from the gain and timing loops. 

Operation of the orthogonal projection operation Pv 1 v 2 x will now be described in relation to the gain and phase 
response of the FIR filter C100. FIG.s 9A and 9B show an example gain and phase response, respectively, for the FIR 
35 filter C100 of the present invention. The gain and phase response vary over time as the filter adapts to parameter 
changes in the recording system; that is, the filter continuously adjusts the channel's overall frequency response so that 
it matches the desired partial response (PR4, EPR4, EEPR4, etc.) as best possible. In the present invention, interfer- 
ence from the timing recovery and gain control loops is attenuated by constraining the gain and phase response of the 
FIR filter C100 at a predetermined frequency using an orthogonal projection operation Pv 1 v 2 _L . 
40 Referring again to FIG. 9A, the gain (magnitude) response of the FIR filter C100 has been constrained to a prede- 
termined value (denoted by g) at the normalized frequency of 0.5 (1/4Ts). Similarly, the phase response of the FIR filter 
C100 has been constrained to kn at the normalized frequency of 0.5 as shown in FIG. 9B. In effect, these constraints 
allow the gain and phase to vary (adapt) at all frequencies except at the normalized frequency of 0.5, thereby constrain- 
ing the filter's frequency response in a manner that attenuates interference from the gain and timing loops. The gain 
45 constraint g is relatively arbitrary except that it is selected to optimize the dynamic range of the filter's coefficients. How- 
ever, constraining the filter's response at the normalized frequency of 0.5 and selecting a phase constraint of kn 
reduces the complexity of the orthogonal projection operation Pv 1 v 2 x > and simplifies implementation of the zero phase 
start circuit B163 (FIG. 4B) and sync mark detector 66 (FIG. 1). 

As mentioned above, the zero phase start circuit B1 63 (FIG. 4B) minimizes the initial phase error between the inter- 
so polated sample values and the baud rate at the beginning of acquisition by computing an initial phase error x from the 
A/D 24 sample values 25 and then loading this initial phase error into the modulo-Ts accumulator B120. To compute 
the initial phase error t, the zero phase start circuit B163 must take into account the phase delay of the adaptive equal- 
izer filter B1 03 since interpolated timing recovery B1 00 operates on the equalized samples 32, not the A/D 24 samples 
25. With the acquisition preamble 68 (FIG. 2B) having a frequency of 1/4T (i.e.. 0.5 normalized), constraining the phase 
55 response of the adaptive equalizer B103 at the preamble frequency (1/4Ts) fixes the phase delay of the equalizer B103 
during acquisition, thereby allowing the zero phase start circuit B163 to accurately compute the initial phase error x. 
Furthermore, since the phase constraint at the preamble frequency is fixed at kn, the phase delay of the equalizer B103 
will either be zero or 180° (i.e., the adjustment to the initial phase error is nothing or a sign change). 

Constraining the phase response of the adaptive equalizer B103 to kn at the preamble frequency also simplifies 
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implementation of the sync mark detector 66 (FIG. 1) in sampled amplitude read channels that use the conventional 
synchronous sampling timing recovery 28 (FIG. 4A) rather than interpolated timing recovery B100 (FIG. 4B). Operation 
of the sync mark detector 66 is described in detail in the above referenced co-pending U.S. Patent Application 
08/533.797 entitled "Improved Fault Tolerant Sync Mark Detector For sampled Amplitude Magnetic Recording." As 
5 described therein, the sync mark detector 66 is enabled coherent with the end of the acquisition preamble 68 and rela- 
tive to the transport delay from the output of the A/D 24 to the sync mark detector 66. With an adaptive equalizer, the 
transport delay will vary unless the filter s phase response is constrained at the acquisition preamble frequency by. for 
example, using an orthogonal projection operation Pv-jVg- 1 of the present invention. 

Turning now to the implementation details of the orthogonal projection operation Pv A v 2 \ the equalizer's frequency 
10 response is 



C(e i2,t, ) = £c, 



-jk2nfT 



20 



where C k are the coefficients of the equalizer's impulse response. At the preamble frequency (1/4T). the equalizer's fre- 
quency response is 



C(e 2 ) = £C k e 



25 where the sampling period has been normalized to T= 1 . In matrix form, the equalizer's frequency response at the pre- 
amble frequency is, 



30 
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40 



.kV 1 



/ n\_(N-l) 







i 


(j)- 1 


= c T 


-j 
-1 


(j) -(N-l) 




j 



Those skilled in the art will recognize that shifting the time base will lead to four different, but functionally equivalent 
frequency responses at the preamble frequency (i.e.. [1.-j,-1.j....]C, [-j.-1 ,j,1....]C, [-1.j.1.-j,..]C and D.Vj.-l .-.]£>. Con- 
45 straining the phase response of the equalizer B103 to an integer multiple of n at the preamble frequency (1/4T) implies 
that the imaginary component of its frequency response is zero, 



so 



55 



0 

-1 

0 

1 



= Q v, = o 



If the imaginary component of the frequency response is constrained to zero, as described above, then constraining the 
magnitude of the equalizer to g at the preamble frequency (1/4T) implies that the real component of the frequency 
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response equals g. 



= C T v 2 = g 



Therefore, the equalizer's coefficients C k must be constrained to satisfy the following two conditions: 

C T -V 1 =0and C T -V 2 = g. 

75 

x The above constraints are achieved by multiplying the computed gradient X k • e^ by an orthogonal projection operation 
Pv 1 v 2 1 as part of a modified LMS algorithm C1 14. 

To understand the operation of the orthogonal projection operation, consider an equalizer that comprises only two 
coefficients: C 0 and C-| as shown in FIG. 10. The phase constraint condition C T * V 1 =0 implies that the filter coeffi- 

20 cient vector C T must be orthogonal to . When using an unmodified LMS algoTithm to update the filter coefficients, the 
orthogonal constraint is not always satisfied as shown in FIG. 10. The present invention, however, constrains the filter 
coefficients to a subspace <C ) which is orthogonal to V-, by multiplying the gradient values X k • e^ by a projection oper- 
ation Pv^, where the null space of the projection operation Pv/ is orthogonal to (£>. The updated coefficients corre- 
spond to a point on the orthogonal subspace <C) closest to the coefficients derived from the unmodified LMS algorithm 

25 as shown in FIG. 10. 

Similar to the phase constraint projection operation Pv^. a second orthogonal projection operation Pv 2 ± constrains 
the filter coefficients such that the coefficient vector C T satisfies the above gain constraint: C T • V 2 = g . The combined 
orthogonal projection operation Pv 1 v 2 - L eliminates two degrees of freedom in an N-dimensional subspace where N is 
the number of filter coefficients (i.e., the orthogonal projection operation Pv 1 v 2 ± has a rank of N-2). 
30 An orthogonal projection operation for V, and V 2 can be computed according to 

Pv/ = I-Pv x = l-V X (V X T V X )- 1 V X T (20) 

where Pv .,v 2 x =Pv } ± • Pv 2 x since V-, is orthogonal to V 2 . The orthogonal projection operation Pv^ 1 computed 
35 using the above equation for an equalizer comprising ten filter coefficients is a matrix 
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SO 



The above matrix Pv 1 v 2 - L is an orthogonal projection matrix scaled by 5 (multiplied by 5) so that it contains integer 
valued elements which simplrf ies multiplying by X k • e^ in the LMS update equation, 

Wm = W k - n - Pv ,v 2 1 • (X k • e k ). (21) 

The scaling factor is taken into account in the selection of the gain value m. Constraining the gain to g and the phase to 
kn at the normalized frequency of 0.5 simplifies implementing the orthogonal projection matrix Pv 1 v 2 1 : half of the ele- 
ments are zero and the other half are either +1.-1, or +4. Thus, multiplying the projection matrix Pv-, v 2 a by the gradient 
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values • ^ requires only shift registers and adders. 

The ACQ/TRK signal shown in FIG. 8B disables adaptation of the FIR filter during acquisition, that is, while acquir- 
ing the acquisition preamble 68 shown in FIG. 2B. Thus, the adaptive equalizer B103 adapts only after acquiring the 
acquisition preamble 68. 

5 

Reduced Cost Orthogonal Constraint Matrix 

Even though the above orthogonal projection matrix Pv, v 2 ± has a simple structure wherein half of the elements are 
zero, it may not be cost effective to directly implement it due to the significant number of shift and accumulate operations 
w necessary to compute Pv 1 v 2 ± • (X k • e k ) . In order to reduce the cost and complexity, an alternative embodiment of 
the present invention decimates the modified LMS adaptation algorithm C1 14 as illustrated in FIG.s 1 1 A 1 1 B 1 1 C and 
11D. 

A mathematical basis for the circuit of FIG. 11 A will be discussed before describing the details of operation. Refer- 
ring again to the above equation (20), 



15 



Pv x X = l-Pv x =l-V x (V x T V x )- 1 V x T 



by combining the above and V 2 vectors into a Nx2 matrix 



20 



25 



V = 



+1 0 
0 -1 
-1 0 
0 +1 



(those skilled in the art will recognize that shifting the time base provides four alternatives for the V matrix) then the 
30 operation (V X T V X )" 1 of equation (20) reduces to 



35 



1 0 
5 1 
0 5 



Thus, equation (20) reduces to 



Pv ± = l-1/5-VV T 



40 



Multiplying both sides of equation (22) by 5 provides 



5*Pv x = 5-l- W T . 

45 Referring again to equation (21), 

W^^W^M-Pv^^-fX.-e,) 
setting X k • e k = g k reduces equation (21) to 

50 

W k+ i =W k -n.Pv x .g k = W H - M .[5-g K -w T .g k ]. 

Defining f> = v T • g k 

55 



(22) 



(23) 



(24) 
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r+i o-io 

[o -1 0 +1 ...J 



go 
Si 

g 3 



8, 



(25) 



70 



15 



SO 



25 



that is, 8 0 =g 0 -g 2 +... and 6 ^-g 1 +g 3 - 
then computing w T • 



+1 0 
0 -1 



-1 0 
0 +1 



and computing 5*g k -w T *g k provides 



30 


+8 0 




5go-5 0 


















35 


+5, 




5g 3 -5 l 



(26) 



(27) 



40 

FIG. 1 1 A implements the above equations (25), (26), (27) and ultimately equation (21) in order to update the coef- 
ficients of the equalizer filter C100 according to the modified LMS adaptation algorithm of the present invention. To fur- 
ther reduce the implementation cost, the circuit of FIG. 1 1 A decimates the adaptation algorithm by the number of filter 
coefficients (10 in the example shown); that is, the adaptation algorithm operates only on every tenth sample value 
45 rather than on every sample value and updates only one filter coefficient per clock period. This is illustrated by expand- 
ing equation (24) 



55 



w 

w (k + 1)0 




W 

w (k + 0)0 




^S(k + 0) 0"^ 


W 

(k + 2) i 




W 




5 g(k + l)l +5 


W 

w (k + 3)2 




W 

w (k + 2)2 


-11- 


5 S(k + 2)2 +5 


w 




w 

W (k + 3)3 




5 S(k + 3)3*^ 



(k+l)l 



(28) 



where k+i is the clock period and g = X H • e s . Thus, at each clock period, the next gradient value g (i)j can be com- 
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puted by multiplying a sample value latched every tenth clock cycle by the next error value ej. The new gradient g, n is 
then used to update the corresponding filter coefficient Wqj of equation (28). 

Referring now to the circuit shown in FIG. 1 1 A, the sample values X k are input into the discrete time equalizer filter 
C100 and the equalized sample values input into the interpolator B122 similar to FIG. 8B. A serial-to-parallel circuit 

5 C140 converts the interpolated sample values Y k+T into even and odd subsequences Y 2k+T and Y^-iw, where the 
notation 2k indicates that two new interpolated sample values are output from the serial-to-parallel circuit 0140 at every 
other sample period. A slicer C142 generates the corresponding even and odd estimated subsequences ~Y 2k+T and 
~ Y (2k-i) + t wh >ch are subtracted from the interpolated sample values at respective adders (0144^01442) to form an 
even and odd sample error sequences e 2k+T and e (2k . 1)+t . An error value interpolation circuit C146, similar to that of 

w FIG. 8B, generates the even and odd sample error sequences e2 k and e (2k -i) which are synchronized to the A/D sample 
rate. 

As mentioned above, the circuit of FIG. 1 1 A decimates the adaptation algorithm of the present invention by the 
number of coefficients in the equalizer filter C1 00. For the example shown, the equalizer filter C100 has 10 filter coeffi- 
cients; accordingly, a decimate by 10 circuit C148 loads a sample value into a delay register C150 every tenth sample 
is period. Thereafter, the output of delay register C150 is represented by X 2k _ 2H where: j=1 ->5 incremented by 2 at every 
other sample period. The output of the delay register C150 X^^, is multiplied by the sample errors e 2k and e (2k „ at 
respective multipliers (C1 52^01 52 2 ) to form the gradient values g 2j and g 2H used in equation (24). 

The gradient values g 2j and g 2H are then shifted into respective shift registers (Cl54 1t C154 2 ). To implement equa- 
tion (27), the gradient values g 2j and g 2H are multiplied by respective alternating ±1 (C156 1f C1562) and accumulated 
so in respective accumulators (C158 1( C158 2 ). After accumulating 5 gradient values in each accumulator (0158^01582), 
the outputs of the accumulators (which represent 5 0 and 8, of equation (25)) are latched by respective decimate by 5 
circuits (0160^01602) and the accumulators (0158^01582) are reset. The values 6q and s, are then multiplied by 
respective alternating +1 (0162^01622) to implement equation (26). The gradient values g 2j and g 2j ^ at the outputs of 
the shift register (C154 1( C154 2 ) are multiplied by 5 (0164^0164^ and the values 6q and ^ are subtracted therefrom 
25 at adders (C166 1 ,C1 66 2 ) in order to implement equation (27). 

To finish the adaptive update algorithm (i.e., to implement equation (28)), the output of adders (C166 1 01 662) are 
scaled by a gain factor u (01 68 t ,01 68 2 ) which is reduced by a factor of 5 to account for the scaled up projection oper- 
ator. The output of the gain factor u (0168^01682) is subtracted at adders (0170^01 70 2 ) from the corresponding filter 
coefficient (W 2j . 1t W 2 0 selected by a multiplexor (0172^01722) from a bank of registers (0174^01742). The ADAPT 2 
signal selects the 2j coefficient from the bank of registers (0174^01742) for updating. After subtracting the update 
value, the updated filter coefficient (C176 1( C176 2 ) is restored to the bank of registers (C174 1( C174 2 ) and used by the 
equalizer filter 01 00 during the next clock period to equalize the sample values according to its updated spectrum 

Decimating the update algorithm by the number of filter coefficients, as in FIG. 1 1 A, decreases the implementation 
complexity but at the cost of slowing convergence of the equalizer filter 01 00 toward an optimum setting. This is 
because the decimated update algorithm does not use all of the sample values to compute the gradients g v FIGs 1 1 B 
and 1 1C show an alternative embodiment of the decimated update algorithm of the present invention wherein more of 
the sample values are used to compute the gradients g j( which adds complexity but improves performance because the 
equalizer filter 0100 converges faster. 

The circuit in FIG. 1 1 B operates similar to the circuit in FIG. 1 1 A described above except for the addition of respec- 
tive gradient averaging circuits (0176^01762) which compute an averaged gradient g 2i , and g 2i over several sample 
values, 1 J 



30 



35 



40 



45 



so 



55 



N- 1 



9 2kj " R 2 e 2K-n * X 2k-2j-n (29) 
n = 0 

where N is a predetermined number of sample values out of the number of equalizer filter taps. In one embodiment, all 
of the sample values could be used (i.e.. a number of sample values equal to the number filter taps), or in an alternative 
embodiment, a decimated number of sample values could be used in order to reduce the complexity and cost. 

FIG. 1 1C shows an embodiment of the gradient averaging circuits (C176 1( C176 2 ) wherein the number of sample 
values is decimated by 2; that is, 5 out of the 10 sample values corresponding to the 1 0 equalizer filter taps are used to 
compute the averaged gradients g 2 j_ t and g 2 j. 

In operation, a decimate by 2 circuit 01 78 of FIG. 1 1 B stores every other sample value in the delay register 01 50 
The output of the delay register 01 50 is multiplied (CISOi-ISOs) by the error value e 2k and delayed versions of the error 
value e^, and the results are accumulated in accumulators (0182^182;;). After accumulating five gradients the con- 
tents of the accumulators (0182^182;;) are transferred to registers (0184^1845), and the accumulators (Cl82 r 182 5 ) 
are cleared. Then, at every other sample period, the contents of registers (Cl84 r 184 5 ) are shifted from left to right (i.e., 
C184 5 = C184 4 ; C184 4 = C184 3 ; etc.) and the output of register C184 5 is the averaged gradient g 2j output by the gra- 
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10 



dient averaging circuit (C176 1( C176 2 ). The averaged gradients g 2j and ^ (0186^01862) are then used to update the 
coefficients of the equalizer filter C100 using circuitry C188 in the same manner as described with reference to FIG 
11A. 

FIG. 1 1 D illustrates yet another embodiment of the present invention which further reduces the cost and complex- 
ity, as compared to FIG. 1 1 A. by updating the even and odd filter coefficients sequentially. That is, during the first N peri- 
ods (where N is the number of fitter coefficients) the circuit of FIG. 1 1 D computes the coefficient update values of 
equation (28) for the even filter coefficients (W 0 , W 2 , W 4 ,..), and then during the next N sample periods it computes the 
update values for the odd filter coefficients (W 1p W 3 , W 5 ,...). A decimation circuit C190 decimates the error value e* by 
two, and approximately half the circuitry as that of FIG. 1 1 A is used to compute the update values. The decimate by two 
circuit C190 is actually syncopated; that is, it outputs the error values for k=0, k=2, k=4, k=6, k=8, and then outputs the 
error values for k=9. k=l 1 , k=13, k=15, k=l7 (assuming the adaptive filter CI 00 comprises ten filter taps). 

Mathematically, operation of the update circuit in FIG. 1 1 D can be described by, 



15 



20 



25 



30 



w 
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W 

w <k+3)2 




W (k + 2)2 


w 

W (k + 5)4 
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w (k + 4)4 


w 

(k+7)6 
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W 

w (k + 9)8 




W 

^ (k + 8)8 


W 

w (k+ 10) 1 




W 

w (k + 9) 1 


W 

w (k+ 12)3 




W 

w (k+ll)3 


W 

W (k+ 14)5 




W (k+ 13)5 


W 

w (k+ 16)7 




W 

w (k + 15)7 


W 

_ (k+ 18)9 




W 

_ (k+ 17)9 



5 g(k + 0)O _5 (ic^0)O 
5 S(k + 2)2 +5 ( k + 2 )0 

5g (k + 4)4 -5 (lc + 4)0 



(k + 6)6 ■ ~(k + 6)0 
>(k + 8)8"^(k + 8)0 
(k + 9)l + 5 ()c + 9)l 
5(k + ll)3" 5 ( k+n >l 
(k + 13)5" h5 ( k + i3)l 



5g 

5g 



(k+l5)7 v (k+15)I 
(k + 17)9 + 5 (k+!7)l 



35 assuming the adaptive filter C1 00 comprises 1 0 filter taps. 

Obviously, the embodiment of FIG. 11 D will decrease the performance of the adaptive filter C100 due to the 
decrease in convergence speed. However, the gradient averaging circuit of FIG. 1 1C can be used to improve the per- 
formance of the or curt in FIG. 1 1D, similar to the embodiment of FIG. 11B. Thus in the embodiments of FIG. 11A-11D, 
the performance versus cost and complexity varies-the preferred configuration is selected according to the require- 

<o ments of the user 

Although the interpolated timing recovery and adaptive equalizer of the present invention have been disclosed in 
relation to a d«0 PR4 read channel, the principles disclosed herein are equally applicable to other types of sampled 
amplitude read channels including d=1 EPR4 or EEPR4 read channels. In a d=1 read channel, for example, the slicer 
B141 of FIG B4A is replaced by a pulse detector as described in the above reference U.S. Patent No. 5,359.631. 
45 Furthermore the particular constraint frequency of 1/4Ts used in the disclosed embodiment is not intended to be 

limiting Other constraint frequencies could be used without departing from the scope of the present invention. For 
example, a 3T preamble could be used in which case the constraint frequency would be 1/6Ts (if constraining to the 
preamble frequency) 

Still further, those skilled in the art will appreciate the many obvious modifications that are possible to the adaptive 
so equations disclosed herein. For example, a shift in the time base would change the FIR filter's magnitude and phase 
response at the constraint frequency which would result in different V 1 and V 2 matrices. Also, the modified LMS algo- 
rithm, 

W k+1 = W h - M .Pv 1 v 2 1 .(X l| .e k ) 

55 

could be implemented after rearranging terms, 

Wk*i = W h . M -e k .(Pv 1 v 2 x .X k ) 
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or 

W k+ i=W k -^.X k -(Pv 1 v 2 ± -e k ). 

5 These, and like modifcations, are within the scope of the present invention. 

The objects of the invention have been fully realized through the embodiments disclosed herein. Those skilled in 
the art will appreciate that the various aspects of the invention can be achieved through different embodiments without 
departing from the essential function. The particular embodiments disclosed are illustrative and not meant to limit the 
scope of the invention as appropriately construed by the following claims. 



Table 1 



Channel 


Transfer Function 


Dipulse Response 


PR4 


(1-DH1+D) 


0,1,0. -1,0, 0, 0, ... 


EPR4 


(1-D) (1+D) 2 


0, 1, 1, -1, -1, 0, 0, ... 


EEPR4 


(1-D)(1+D) 3 


0,1,2, 0,-2,-1,0, ... 
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Table B2 



5 



10 



25 



35 



40 



50 



x-32/Ts 


c(-2) 


C(-2) 


C(0) 


C(1) 


c(2) 


c(3) 


0 


0.0000 


-0.0000 


1.0000 


0.0000 


-0.0000 


0.0000 


1 


0.0090 


-0.0231 


0.9965 


0.0337 


-0.0120 


0.0068 


2 


0.0176 


-0.0445 


0.9901 


0.0690 


-0.0241 


0.0135 


3 


0.0258 


-0.0641 


0.9808 


0.1058 


-0.0364 


0.0202 


4 


0.0335 


-0.0819 


0.9686 


0.1438 


-0.0487 


0.0268 


5 


0.0407 


-0.0979 


0.9536 


0.1829 


-0.0608 


0.0331 


6 


0.0473 


-0.1120 


0.9359 


0.2230 


-0.0728 


0.0393 


7 


0.0533 


-0.1243 


0.9155 


0.2638 


-0.0844 


0.0451 


8 


0.0587 


-0.1348 


0.8926 


0.3052 


-0.0957 


0.0506 


9 


0.0634 


-0.1434 


0.8674 


0.3471 


-0.1063 


0.0556 


10 


0.0674 


-0.1503 


0.8398 


0.3891 


-0.1164 


0.0603 


11 


0.0707 


-0.1555 


0.8101 


0.4311 


-0.1257 


0.0644 


12 


0.0732 


-0.1589 


0.7784 


0.4730 


-0.1341 


0.0680 


13 


0.0751 


-0.1608 


0.7448 


0.5145 


-0.1415 


0.0710 


14 


0.0761 


-0.1611 


0.7096 


0.5554 


-0.1480 


0.0734 


15 


0.0765 


-0.1598 


0.6728 


0.5956 


-0.1532 


0.0751 


16 


0.0761 


-0.1572 


0.6348 


0.6348 


-0.1572 


0.0761 


17 


0.0751 


-0.1532 


0.5956 


0.6728 


-0.1598 


0.0765 


18 


0.0734 


-0.1480 


0.5554 


0.7096 


-0.1611 


0.0761 


19 


0:0710 


-0.1415 


0.5145 


0.7448 


-0.1608 


0.0751 


20 


0.0680 


-0.1341 


0.4730 


0.7784 


-0.1589 


0.0732 


21 


0.0644 


-0.1257 


0.4311 


0.8101 


-0.1555 


0.0707 


22 


0.0603 


-0.1164 


0.3891 


0.8398 


-0.1503 


0.0674 


23 


0.0556 


-0.1063 


0.3471 


0.8674 


-0.1434 


0.0634 


24 


0.0506 


-0.0957 


0.3052 


0.8926 


-0.1348 


0.0587 


25 


0.0451 


-0.0844 


0.2638 


0.9155 


-0.1243 


0.0533 


26 


0.0393 


-0.0728 


0.2230 


0.9359 


-0.1120 


0.0473 


27 


0.0331 


-0.0608 


0.1829 


0.9536 


-0.0979 


0.0407 


28 


0.0268 


-0.0487 


0.1438 


0.9686 


-0.0819 


0.0335 


29 


0.0202 


-0.0364 


0.1058 


0.9808 


-0.0641 


0.0258 


30 


0.0135 


-0.0241 


0.0690 


0.9901 


-0.0445 


0.0176 


31 


0.0068 


-0.0120 


0O337 


0.9965 


-0.0231 


0.0090 



55 Claims 



A sampled amplitude read channel for reading digital data from a sequence of discrete time sample values gener- 
ated by sampling an analog read signal from a read head positioned over a magnetic medium, comprising: 
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(a) a sampling device for sampling the analog read signal to generate the discrete time sample values; 

(b) an adaptive equalizer, responsive to the discrete time sample values, for generating equalized sample val- 
ues according to a target response; 

(c) interpolated timing recovery for generating interpolated sample values; and 

5 (d) a discrete time sequence detector for detecting the digital data from the interpolated sample values. 

2. The sampled amplitude read channel as recited in claim 1, wherein: 

(a) the adaptive equalizer comprises an FIR filter comprising a plurality of filter coefficients; 

10 (b) the adaptive equalizer operates according to a least mean square algorithm responsive to an error value ^ 

computed as a function of an output of the FIR filter and an estimated ideal value; 

(c) the adaptive equalizer comprises an interpolation circuit, responsive to the interpolated sample values, for 
generating the error value synchronous with the discrete time sample values; and 

(d) the least mean square algorithm uses the error value e k to adjust the filter coefficients of the FIR filter. 

15 

3. The sampled anplitude read channel as recited in claim 2, wherein the error value is computed relative to a dif- 
ference between an output of the interpolated timing recovery and a sample value estimator responsive to the inter- 
polated sample values. 

20 4. The sampled amplitude read channel as recited in claim 2, wherein the error value ^ is computed relative to a dif- 
ference between an output of the FIR filter and an output of the interpolation circuit 

5. The sampled amplitude read channel as recited in claim 2, wherein: 
25 (a) the least mean square algorithm is: 

(b) Wj. are the filter coefficients of the FIR fSter; 
30 (c) 1 1 is a predetermined gain; and 

(d) ts a discrete time sample value. 

6. The sampled amplitude read channel as recited in claim 2, wherein the adaptive equalizer further comprises a con- 
straint circuit tor constraining a frequency response of the adaptive equalizer at a predetermined constraint fre- 

35 quency 

7. The sampled amplitude read channel as recited in claim 6, wherein the constraint circuit comprises an orthogonal 
projection operation tor constraining the filter coefficients relative to the predetermined constraint frequency. 

40 8. The sampled amplitude read channel as recited in claim 7, wherein: 

(a) the constraint circuit operates according to a modified least mean square algorithm, 

W k+ i = W k - M -Pv 1 v 2 - L -(X k .e k ) 

45 

(b) the filter coefficients; 

(c) |i is a predetermined gain; 

(d) Pv,vj is the orthogonal projection matrix; and 

(e) X* is a discrete time sample value. 

so 

9. The sampled amplitude read channel as recited in claim 1, wherein the interpolated timing recovery further gener- 
ates a data clocfc tor clocking the discrete time detector. 

10. The sampled amplitude read channel as recited in claim 1 , wherein: 

55 

(a) the analog read signal comprises pulses modulated by the digital data at a predetermined baud rate; and 

(b) the sampling device samples the analog read signal at a rate above the baud rate. 

11. The sampled amplitude read channel as recited in claim 1 , wherein the interpolated timing recovery comprises an 
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interpolator responsive to an interpolation interval % proportional to a time difference between a channel sample 
value and a desired interpolated sample value. 

12. The sampled amplitude read channel as recited in claim 1 , wherein the analog read signal comprises pulses mod- 
5 ulated by the digital data at a predetermined baud rate, the interpolated timing recovery further comprising; 

(a) a frequency offset generator for generating a frequency offset Af proportional to a frequency difference 
between a sampling clock and the baud rate; and 

(b) a mod-Ts accumulator for accumulating, modulo-Ts, the frequency offset Af to generate an interpolation 
io interval t where Ts is a predetermined sample period of the sampling clock. 

13. The sampled amplitude read channel as recited in claim 12, wherein the frequency offset generator comprises: 

(a) a phase error detector for detecting a phase error AO between an interpolated sample value and an esti- 
is mated sample value; and 

(b) a loop filter for filtering the phase error AO to generate the frequency offset Af, 

14. The sampled amplitude read channel as recited in claim 11, wherein the interpolator comprises an actual interpo- 
lation filter having an actual impulse response h(k) which approximates an ideal impulse response of an ideal inter- 
ne polation filter: 

sincfrc * (k-x/Ts)); where: 

k is a time index; and 
25 Ts is a sample period of the sampling clock. 

15. The sampled amplitude read channel as recited in claim 11, wherein the interpolator comprises an actual interpo- 
lation filter having an actual impulse response h(k) generated by minimizing a mean squared error between a fre- 
quency response of the actual interpolation filter and an ideal frequency response corresponding to an ideal 

30 interpolation filter. 

16. The sampled amplitude read channel as recited in claim 1 1 , wherein: 

(a) the interpolator comprises an actual interpolation filter; and 

35 (b) a plurality of coefficients C T k of the actual interpolation filter are computed in real time as a function of t. 

17. A method for reading digital data from a sequence of discrete time sample values generated by sampling and ana- 
log read signal from a read head positioned over a magnetic medium, comprising: 

4 o (a) sampling the analog read signal to generate the discrete time sample values; 

(b) adaptively equalizing the discrete time sample values to generate equalized sample values according to a 
target response; 

(c) interpolating the equalized sample values to generate interpolated sample values; and 

(d) detecting the digital data from the interpolated sample values. 

45 

18. The method for reading digital data as recited in claim 17, wherein the step of interpolating sample values com- 
prises the step of generating a data clock for clocking the discrete time sequence detector. 
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